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(57) Abstract 

As stereotactic images are digitised, the method consists in the following steps: 
selecting (400) a target pixel in a target region of interest, generating (401), around the 
selected pixel, a target window of selected dimensional characteristics and containing the 
said target region of interest, determining (402) in a second image, based on a prcdeicrmincd 
criterion of selection, a set of pixels, generating (403) around each selected pixel a second 
window of the same dimensional characteristics as the said target window, effecting a 
correlation processing (404) between the grey levels of the pixels of each second window 
and the grey levels of the target window pixels, so as to obtain a correlation value for each 
second window, and identifying the region of interest homologous with the target region 
of interest based on the analysis of the whole set of coirelation values obtained, so as to 
minimise the risks of pairing errors between homologous regions of interest. Then the 
element of interest is located from the positions of the two homologous regions. 

(57) Abr^g^ 

Les images st6r6otaxiques 6tant num6ris6es, on s61ectionne (400) un pixel cible dans 
une r6gion d'int^ret cible, on genfcre (401), autour du pixel cible s61ectionn€, une fenetre 
cible de caractdristiques dimensioiuielles choisies et contenant ladite region d'int^r^t cible, 
on ddtermine (402) dans imc deuxifeme image, selon un entire de selection predetermine, 
un ensemble de pixeb, on gdn^re (403) autour de chaque pixel sdlectionne une deuxi^me 
fenfitre de m&mes caractdristiques dimensionnelles que ladite fen6tre cible, on effectue un 
traitement de conflation (404) entre les niveaux de gris des pixels de chaque deuxi&me 
fenetre et les niveaux de gris des pixels de la fenfttre cible, de fa^on ^ obtenir une valeur de 
correlation poiir chaque deuxi&me fenetre, et on identifie la region d'inter€t homologue de 
la region d*inter&t cible ^ partir de I'analyse de I'ensemble des valeurs de correlation ainsi 
obtenues, de fa^on 2h minimiser les risques d'erreur d'appariement entre les regions d'intergt 
homologues. On localise ensuite reiernent d'inter&t ^ partir des positions des deux regions homologues. 
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Procede de localisation d'un element d'interet contenu dans un 
objet tridimensionnel, en particulier lors d*un examen de 
stereotaxic en mammographie. 

L'invention conceme la localisation d'un element d'interet contenu 
dans un objet tridimensionnel a partir des positions de regions 
d'interet homologues correspondant audit element d'interet et 
apparaissant dans un ensemble d'images stereotaxiques dudit objet. 

Bien que I'invention s'applique h tout ensemble d'images 
stereotaxiques d'un objet tridimensionnel quelconque, elle est 
particulierement utile dans le domaine medical et notamment lors d'un 
examen de stereotaxic en mammographie pour la localisation de 
microcalcifications dans un sein. 

Un examen de stereotaxic en general jpermet d'acc6der avec une 
grande precision, de Tordre du millimetre, h un point donne d'un objet 
tridimensionnel, a partir de deux projections bidimensionnelles de cet 
objet acquises par exemplc selon deux incidences angulaires opposees 
de part et d'autre de la normale au plan du receptcur d'images delivrant 
lesdites projections. 

A partir de la connaissance, d'une part, de la position de la 
projection de ce point dans le plan dc chaque image stereotaxique 
obtenue, et, d'autre part, de la geometric de I'appareil de prises de 
vues stereotaxiques ayant conduit a la constmction des deux images, il 
est possible de calculer la position exacte spatiale de ce point dans 
Tobjet tridimensionnel, par calcul trigonom6trique. 

En resume, un examen stereotaxique necessite : 

de faire au moins deux images de Tobjet k differentes angulations, 
de connaitre parfaitement la geom6trie d'acquisition du systeme 
stereotaxique, et 

de pouvoir localiser sur les differentes images stereotaxiques 
obtenues, le lieu de la projection de r61ement d'interet choisi. 

Les deux premiers points ne posent pas de problfeme majeur. 

Actuellement, en mammographie, rexamen de stereotaxic est 
effectu6 k Taide d'un mammographe equip6 d'un dispositif de prises de 
vues stereotaxiques, Le manrunographe comprend un tube k rayons X, 
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situe a rextremite d'un premier bras qui est mobile autour d'un axe et 
emettant un rayonnement X vers un recepteur situe a Textremite d'un 
autre bras, Entre le tube et le recepteur sont disposes, d'une part, un 
plateau porte-sein, ou support patient, et, d'autre part, un plateau de 
compression assurant le maintien en place du sein lors des 
mammographies. Le recepteur d'images peut 6tre un recepteur 
numerique, tel qu'une camera CCD par exemple, delivrant des images 
stereotaxiques numerisees. Les prises de vues stereotaxiques 
necessitent la rotation du tube a rayons X autour des plateaux porte- 
sein et de compression, selon deux orientations successives opposees 
de part et d'autre de sa position initiale perpendiculaire au plan du 
recepteur d'images. 

A partir des images stereotaxiques numerisees obtenues, la 
localisation de I'element dlnterSt sur les differentes images est une 
operation delicate, en particulier en mammographie, puisqu'elle fait 
appel a une nndse en correspondance, ou appariement ("matching" en 
langue anglaise) de regions d*interet homologues correspondant, dans 
deux images stereotaxiques differentes, a un meme element d'interet 
situe dans le sein. 

Ainsi, en mammographie, on cherche a apparier sur les differentes 
images stereotaxiques obtenues, les regions d'interet qui correspondent 
a une meme microcalcification. 

Or, I'appariement des microcalcifications projetees est un 
probleme difficile car les microcalcifications ne se ressemblent pas 
d'une image a Tautre. Leur forme et leur contraste peuvent changer 
ainsi que leur agencement dans I'espace. Elles peuvent se superposer 
avec des zones fibreuses ou avec d*autres microcalcifications. 

Jusqu'a present, Tappariement etait effectue entierement 
manuellement par le radiologue. Ce demier choisissait quelles etaient 
les regions d'interSt en correspondance sur les deux images 
stereotaxiques angulees, puis rentrait les coordonnees 
bidimensionnelles de ces regions dans un calculateur qui donnait la 
position spatiale de la microcalcification, c*est-a-dire les coordonnees 
tridimensionnelles de cette demiere. A partir de ces coordonnees, le 
radiologue pouvait ajuster par exemple la position d*un porte-aiguille 
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de fagon a ce que I'aiguille vienne a Taplomb de la microcalcification 
pour permettre sa ponction par exemple. 

Or, cet appariemment entierement manuel s'avere etre une 6tape 
longue et parfois imprecise. 
5 L'invention vise k apporter une solution k ce probleme. 

L'invention a pour but de reduire sensiblement la duree d'un 
examen stereotaxique et d'obtenir une meilleure precision 
d'appariement. 

D*une fa9on tres generate, selon Tinvention, on definit une mesure 

10 de similarite permettant d'evaluer la ressemblance entre une region 

d'interet (microcalcification par exemple) selectionnee (manuellement 
par le radiologue par exemple, ou bien automatiquement par un 
procede de selection automatique) sur une des images stereotaxiques 
angulees et une region d'interet d'une deuxieme image stereotaxique. 

15 A partir de cette mesure de similarite, on selectionne, pour obtenir 

la region d'interet homologue, un certain nombre de "candidats" 
appartenant a une zone de recherche et ayant une forte ressemblance 
avec la microcalcification de depart. 

Ensuite, en utilisant les principes de la stereotaxic, on calcule les 

20 coordonnees tridimensionnelles de ces candidats, puis on les reprojette 

sur I'image de centrage. On calcule a nouveau la similarite entre ces 
points reprojetes et la microcalcification initialement choisie. En 
effet, si le candidat correspond effectivement a la microcalcification 
de depart, alors sa projection calculee correspondra effectivement a 

25 rimage de cette naicrocalcification sur Fimage de centrage. 

Enfin, a Taide de ces deux valeurs de similarite calculees pour 
chaque candidat, on propose le meilleur candidat. 

En theorie, on pourrait s'abstenir de I'etape de reprojection 
utilisant I'image de centrage et se contenter de Tappariement tel qu'il 

30 vient d'etre expos6 pour selectionner la region d'interet homologue de 

celle selectionn6e sur I'image cible. Cependant, il est nettement 
preferable d'utiliser la reprojection des candidats sur I'image de 
centrage pour verifier la premiere s61ection effectuee, notamment dans 
le domaine de la msimmographie ou les microcalcifications peuvent 

35 changer de forme et de contraste d'une image a I'autre. 
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En d'autres termes, Tinvention propose tout d'abord un proc6de de 
localisation d'un 61ement d'int6ret contenu dans un objet 
tridinciensionnel a partir des positions de regions d'interSts homologues 
correspondant audit Element d*interet et apparaissant dans un ensemble 
d'images stereotaxiques dudit objet, Ce procede conciporte une etape de 
selection dans une premiere image ster^otaxique, ou image cible, d'une 
premiere region d'interSt, en particulier une microcalcification, dite 
region d'interet cible. Le procede comporte egalement un appariement 
de ladite premiere region avec une deuxieme region d'interet 
homologue de la premiere et apparaissant dans une deuxieme image 
stereotaxique. 

Selon une caracteristique generale de I'invention, les images 
stereotaxiques etant numerisees, on selectionne dans I'etape de 
selection un pixel cible dans ladite region d'interet cible. Dans Tetape 
d'appariement, on genere autour du pixel cible selectionne, une fenetre 
cible de caracteristiques dimensionnelles choisies et contenant ladite 
region d'int6r&t cible. On determine alors dans la deuxieme image, 
selon un critere de selection pr6determine, un ensemble de pixels et on 
genere autour de chaque pixel selectionne une deuxieme fenetre de 
memes caracteristiques dimensionnelles que ladite fenetre cible. On 
effectue un traitement de correlation entre les niveaux de gris des 
pixels de chaque deuxieme fenetre et les niveaux de gris des pixels de 
la fenetre cible de fa9on a obtenir une valeur de correlation pour 
chaque deuxieme fenStre. On identifie alors la region d'interet 
homologue de la region d'interSt cible a partir de Tanalyse de 
Tensemble des valeurs de correlation ainsi obtenues. On minimise 
ainsi les risques d'erreur d'appariement entre des regions d'interSt 
homologues. 

Selon un mode de mise en oeuvre, Tanalyse des valeurs de 
correlation obtenues comporte la selection d'un nombre choisi de 
maxima ou de minima de correlation, la selection de la region 
homologue s'effectuant parmi celles dont la valeur de correlation 
associee est Tun de ces maxima ou minima de correlation. 

En fait suivant que I'image est "normale" ou inversee (une image 
"normale" etant definie comme presentant des niveaux de gris et des 
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contrastes qui sont semblables a ceux habituellement visibles sur des 
cliches de film) on utilisera respectivement les maxima ou les minima. 

On peut egalement prevoir une detemiination de la dynamique des 
maxima ou minima selectionnes, et la comparaison des valeurs de 
dynamique obtenues k un seuil. 

Le traitement de correlation utilis6 peut comporter un traitement 
de correlation normee, centree ou non, ou bien un traitement de 
difference normee. 

En variante a I'appariement evoque ci-dessus, Tinvention propose 
egalement un appariement base sur Telaboration de caracteristiques 
pr6determinees de toute region d'interet, en particulier des 
caract6ristiques physiques comme Felongation, le contraste, le 
gradient, puis une comparaison entre ces listes de caracteristiques 
associees a differentes regions dlnterets. 

Plus precisement, selon une caracteristique generale de 
rinvention, on s61ectionne dans T^tape de selection, un pixel cible 
dans ladite region d*interet cible et, dans Tetape d'appariement, on 
genere, autour du pixel cible selectionne, une fenetre cible de 
caracteristiques dimensionnelles choisies et contenant ladite region 
d'interet cible. On determine pour ladite region d'interSt cible, a partir 
des valeurs de niveau de gris des pixels de ladite fenetre cible, un 
premier ensemble de valeurs numeriques respectivement 
representatives de caracteristiques predeterminees caracterisant toute 
region d'interet (par exemple des caracteristiques de forme de 
contraste, de gradient ... ) et on stocke ce premier ensemble de valeurs 
numeriques. On determine dans la deuxieme image, selon un critere de 
selection predetermine, un ensemble de pixels appartenant 
respectivement a des regions dlnterets dites candidates apparaissant 
dans la deuxieme image et on genere autour de chacun de ces pixels 
selectionnes une deuxieme fenStre contenant ladite region d*interet 
candidate correspondante. On detemndne, pour chaque region d'interet 
candidate a partir des valeurs de niveau de gris des pixels de la 
deuxieme fenetre associee, un deuxieme ensemble de valeurs 
numeriques respectivement representatives desdites caracteristiques 
pr6d6termin6es et on stocke ce deuxifeme ensemble de valeurs 
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numeriques. On identifie alors la region d^interet homologue de la 
region d'int6ret cible a partir d'un traitement de comparaison entre le 
premier ensemble de valeurs numeriques et chaque deuxieme ensemble 
de valeurs numeriques. On minimise egalement ainsi les risques 
5 d'erreur d'appariement entre des regions homologues. 

Le traitement de comparaison peut comporter un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

Dans les deux variantes d'appariement qui viennent d'etre 
10 evoquees, et afin de nciinimiser encore le temps de calcul, on determine 

dans la deuxieme image stereotaxique, une zone epipolaire contenant 
au moins le segment epipolaire relatif au pixel cible. Les pixels de la 
zone epipolaire contiennent alors les pixels selectionnes. Pour reduire 
encore le temps de calcul, au lieu d'effectuer Tappariement sur tous les 
15 pixels de la zone epipolaire determinee, on selectionne 

avantageusement les pixels de la deuxieme image parmi ceux 
presentant des maxima ou des minima de niveau de gris dont la 
dynamique est superieure a un seuil predetermine. 

En outre, pour s'affranchir des problemes de fond d'image, et pour 
20 minimiser encore les risques d'erreurs d'appariement entre deux 

regions dlnterets homologues, il est particulierement avantageux que 
Tappariement selon Tinvention comporte prealablement un filtrage des 
images stereotaxiques, par exemple du type dit "chapeau haut de 
forme". 

25 L'utilisation de Tetape de reprojection des elements d*interets 

candidats dans Timage de centrage est independante du type 
d'appariement utilise prealablement. 

Ainsi, selon une caracteristique generale de Tinvention, le procede 
de localisation du type de celui evoque ci-avant comprend une 

30 selection dans une premiere image stereotaxique numerisee d*une 

premiere region d'interet dite cible, une selection, dans une deuxieme 
image stereotaxique numerisee, a partir d'un prennder appariement 
automatique, d'au moins une deuxieme region d'interet susceptible 
d'etre homologue de la region cible, la determination de la position 

35 spatiale d'un pixel candidat d'un element d'interet candidat 
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correspondant a ces deux regions d'interets, la d6termination de la 
position dans une troisieme image stereotaxique, d'un pixel projete 
correspondant a la projection dans cette troisieme image stereotaxique 
du pixel candidat, et un deuxieme appariement automatique entre la 
region interet cible et un voisinage du pixel projete definissant une 
region d'intergt projete. 

En pratique, on selectionne avantageusement dans la deuxieme 
image stereotaxique, a partir dudit premier appariement automatique, 
un ensemble de deuxiemes regions susceptibles d'etre homologues de 
la region cible et on determine la position spatiale de chaque pixel 
candidat correspondant a chaque paire de regions d'interSts formSe par 
la region cible et Tune des deuxiemes regions. On determine 
egalement la position dans la troisieme image de chaque pixel projete 
correspondant et on effectue le deuxieme appariement automatique 
entre la region cible et un voisinage de chaque pixel projete. 

On pent ordonner Tensemble des deuxiemes regions obtenues dans 
I'ordre decroissant de leur probabilite a etre la region homologue de la 
region cible. On peut egalement ordonner Tensemble des regions 
projetees obtenues dans Tordre decroissant de leur probabilite h Stre la 
region homologue de la region cible. On selectionne alors 
avantageusement la region homologue de la region cible comme etant 
celle dont le produit des rangs respectifs dans les deux 
ordonnancements est mininum. 

Le prennier appariement automatique peut etre Tun des deux 
appariements qui ont ete evoques ci-avant, c'est-a-dire a base d'un 
traitement de correlation ou d'un traitement de comparaison. 

De meme, le deuxieme appariement automatique peut comporter 
un traitement de correlation tel que defini ci-avant, entre le voisinage 
de chaque pixel projete et une fenStre cible contenant la region 
d'interet cible. Ce deuxieme appariement automatique peut 6tre 
egalement un traitement de comparaison tel que d6fini ci-avant, entre 
la region d'interet cible et chaque region d'int6ret projetee. 

D'un point de vue materiel, Tinvention s'analyse comme un 
dispositif de localisation d'un element d'interet contenu dans un objet 
tridimensionnel a partir des positions de regions d'interets homologues 
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correspondant audit element d*int6ret et apparaissant dans un ensemble 
damages stereo taxiques dudit objet. Ce dispositif comporte un 
appareillage de prises de vues stereotaxiques equipe d'un support pour 
ledit objet et d*un recepteur damage tel qu'une camera CCD capable de 
foumir des images ster6otaxiques numeris6es. II est egalement prevu 
des moyens de visualisation de ces images, par exemple un ecran 
video, tel que celui d'un microordinateur. Des moyens de selection, 
par exemple une "souris" d*un microordinateur permettant de deplacer 
un marqueur sur TScran et de "cliquer" ce marqueur sur un pixel cible 
d'une region d'interet, permettent la selection du pixel cible. II est 
enfin prevu un microprocesseur et des memoires associees' pour 
incorporer sous forme logicielle des moyens permettant de realiser 
fonctionnellement toutes les etapes des procedes qui ont 6t6 decrites 
ci-avant. 

D*autres avantages et caracteristiques de Tinvention apparaitront a 
Texamen de modes de mise en oeuvre nullement limitatifs, et des 
dessins annexes, sur lesquels : 

la figure 1 est une vue schematique de trois images stereotaxiques 
d'un objet tridimensionnel a partir d'un dispositif de prises de vues 
stereotaxiques ; 

la figure 2 illustre tres schematiquement la notion de segment 
epipolaire; 

la figure 3 illustre une zone epipolaire choisie conformement a 
Tinvention pour effectuer un appariement; 

la figure 4 est un organigramme general d'un procede selon 
rinvention; 

la figure 5 illustre plus particulierement une 6tape d'une variante 
dii procede utilisant un traitement de correlation; et 

la figure 6 illustre schematiquement plus particulierement un 
procede selon Tinvention utilisant une reprojection sur Timage de 
centrage. 

Un examen stereotaxique, en particulier en mammographie, se 
compose d'une serie de trois expositions d'un objet tridimensionnel 1 
(figure 1), par exemple un sein, reposant sur un support 3 et comprime 
par un plateau de compression 2, a I'aide d'un tube a rayons X 
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occupant respectivement trois positions differentes 6, 7 et 8. En 
pratique, on procede k une exposition selon un angle de 0° et a deux 
expositions angulees selon deux angles opposes egaux, en pratique 
±15°. 

Un recepteur CCD 4 permet d*obtenir ainsi trois images 
stereotaxiques numerisees 10, 11 et 12. Selon une convention 
habituellement utilisee, Timage 10 est Timage de droite, tandis que 
I'image 12 est Timage de gauche et que Timage 11 est Timage de 
centrage. 

Un element d*interet 5, par exemple une microcalcification, 
contenu dans Tobjet tridimensionnel 1, foumit sur chacune des ifnages 
10, 11 et 12 des regions d'interets respectivement referencees 5d, 5c et 
5g. L'image de centrage sert en particulier a verifier que la 
microcalcification sera accessible durant Facte medical eventuel 
ulterieur. 

Les coordonnees tridimensionnelles de la microcalcification 5 
peuvent Stre obtenues, de fa9on classique et bien connue de I'homme 
du metier, par un calcul simple de trigonometric a partir de la 
connaissance des coordonnees bidimensionnelles des deux projections 
5d et 5g, par exemple. Ces coordonnees sont exprimees par rapport a 
une origine materialisee par le centre d'une croix 9 gravee sur le 
plateau supportant le sein. Et on repere sur chaque image le centre de 
cette croix de reference pour initialiser geometriquement le systeme. 

L'homme du metier sait que, par definition, le segment epipolaire 
50 (figure 2) d*un point 5d d'une image 9 d'angulation 61, sur une 
image 12 d'angulation 62, correspond a Tensemble des points de la 
seconde image 12 pouvant avoir le point 5d comme antecedent. 
Aiatrement dit, c'est la projection conique par la source d'angulation 
92, du segment de droite 52 joignant ledit point a la source 
d'angulation 61. 

Dans le cas oii les images stereotaxiques ont etc obtenues avec des 
angulations telles que les positions du tube a rayons X se situent dans 
un plan parallele aux lignes des images, ce segment epipolaire 50 est 
un segment de droite qui se trouve etre parallele aux lignes de Timage. 
Cependant, si Ton n'a pas la connaissance de I'epaisseur de I'objet 
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tridimensionnel, on ne peut pas determiner la position des extremites 
de ce segment epipolaire. Ce segment epipolaire devient alors une 
droite epipolaire qui, dans la geom6trie d'acquisition particuli&re 
evoquee ici, se situe sur la meme ligne d'image 12 que le point 

5 d'interet 5d sur Timage 9. 

En outre, lorsqu*on utilise un detecteur qui presente des 
distorsions geometriques, la droite epipolaire se transforme en une 
courbe epipolaire 53 (figure 3). On peut determiner les coordonnees de 
I'ensemble des points de la courbe epipolaire sur Timage 12 par une 

10 calibration geometrique de la camera CCD. II serait alors possible de 

n*effectuer Tappariement selon I'invention que sur I'ensemble de ces 
points. 

Cependant, si pour des raisons de precision de la calibration 
geometrique, on ne peut determiner qu'approximativement les 
15 coordonnees des points de la courbe epipolaire, on choisira alors 

d'effectuer Tappariement selon Tinvention sur un ruban epipolaire 51 
qui encadre la courbe epipolaire. En pratique, on definira une zone 
epipolaire, ou zone de recherche, delimitee par deux lignes parall&les 
passant respectivement par deux points augmentes de part et d'autre du 
20 point le plus haiit et le plus bas de Tepipolaire 53 d'une marge 

predeterminee. Generalement, on opte pour une marge de 3 a 5 pixels. 

On va maintenant decrire plus en detail, en se referant plus 
particulierement aux figures 4 et 5, deux types d'appariement selon 
rinvention. 

25 La premiere etape consiste en la selection (etape 400) d'un pixel 

cible 60d au sein d'une region d'interSt cible 5d presente par exemple 
dans rimage de droite 10, dite image cible. En pratique, cette 
selection peut s'effectuer manuellement en cliquant un point de la 
region d*interSt de I'image concemee apparaissant sur un ecran d*un 

30 microordinateur. Ceci etant, on peut prevoir egalement un precede de 

selection automatique de pixel. Ce pixel 60d peut etre centre ou non 
dans la region cible. 

On genere ensuite (etape 401) une fenetre cible 61 d contenant 
ladite r6gion d'interet cible. Cette fenStre cible peut etre gen^ree 

35 automatiquement apres avoir applique sur le pixel cible 60d un moyen 
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logiciel classique de croissance de region. En variante,^ on peut 
egalement prevoir, selon le domaine d'application de Tinvention, de 
choisir des fenetres de taille pred6finies. Dans tous les cas, on 
obtiendra generalement une fenetre cible rectangulaire ou carree dont 

5 les c6tes presentent des dimensions de Tordre d'une dizaine a quelques 

dizaines de pixels. 

Les niveaux de gris de chacun des pixels de cette fenetre cible 
sont ensuite stockes dans une memoire associee a un microprocesseur 
qui incorpore de fa9on logicielle tous les moyens fonctionnels de mise 

10 en oeuvre du procede decrit. 

A partir de Ik, on va selectionner dans une etape 402 un ensemble 
de pixels dans une deuxi^me image stereotaxique, par exemple Timage 
de gauche. Get ensemble de pixels selectionnes va constituer une zone 
de recherche 51 a Tinterieur de laquelle on va selectionner une ou 

15 plusieurs regions dlnterets candidates susceptibles d'etre homologues 

de la region d'interet cible. 

Comme on Ta explique ci-avant et dans le but de reduire la duree 
d'appariement, cette zone de recherche 51 sera par exemple la zone 
epipolaire 51 definie ci-avant. Plus precisement, autour de chacun de 

20 ces pixels selectionnes dans la zone de recherche, on genere (etape 

403) une deuxieme fen6tre 61g de meme caracteristique 
dimensionnelle que la fenetre cible. On stocke les niveaux de gris des 
pixels de chaque deuxieme fenStre et on procede (etape 404) a un 
traitement de correlation entre les niveaux de gris des pixels de chaque 

25 deuxieme fenetre 61 g et les niveaux de gris des pixels de la fenetre 

cible 6 Id. 

Le traitement de correlation permet done d*effectuer une mesure 
de similarite entre la fen6tre cible et chacune des deuxiemes fenetres 
de la deuxieme image de fa9on k obtenir une valeur de similarit6 
30 (valeur de correlation entre ces deux fenStres). 

C'est a partir de I'analyse de ces valeurs de correlation que Ton va 
proceder a I'appariement de la r6gion d'interet cible avec sa region 
d'interSt homologue. 

En ce qui conceme le traitement de correlation (6tape 404), on 
peut utiliser une correlation norm6e. Plus pr6cisement, chaque valeur 
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de correlation d*une deuxieme fenStre sera alors donnee par la 
formule : 

N 

^ I (k) J (k) ^ 
k=i I.J 



V^i J (k) mil II J II 

ou N est le nombre de pixels d'une fenetre, k Tindice courant d*un 
10 pixel, I(k) le niveau de gris du pixel k de la fenetre cible et J(k) le 

niveau de gris du pixel k de la deuxieme fenetre. 

On pent egalement utiliser une difference normee permettant 
d'obtenir la valeur de correlation par la formule : 



II T 



En theorie, la region d'interet homologue sera celle contenue dans 
la deuxieme fenetre qui presente la valeur de correlation maximale. 
Ceci etant, en pratique, plusieurs maxima de correlation peuvent etre 

20 trouves avec une intensite elevee. On selectionne alors (etape 405), 

parmi I'ensemble des valeurs de correlation, un ensemble restreint de 
maxima de correlation selon un critere predetermine comme Tintensite 
du maximum ou comme le nombre de maxima dans la liste des valeurs 
de correlation rangees par ordre dlntensite d6croissante. 

25 On peut egalement utiliser coxrune critere supplementaire de 

selection une d6termination de la dynamique des maxima selectionnes 
et la comparaison des valeurs de dynamique obtenues a un seuil 
predetermine. 

La determination de la dynamique de maxima est parfaitement 
30 connue de Thomme du metier. Cependant, celui-ci pourra 

eventuellement se reporter a la demande de brevet fran9ais n° 
91 15308 pour plus de d6tails concemant une telle d6termination de 
dynamique, 

II convient de noter ici que le procede selon I'invention permet de 
35 determiner et d'afficher sur r6cran du nriicroordinateur, par exemple 
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par un clignotement, la region d'interet consid6r6e comme etant la 
region homologue de la region cible, et ce avec le minimum de risque 
d'erreur. 

Ceci etant, en pratique, le procede selon I'invention pourra 

5 selectionner un nombre restreint de regions d'interSts candidates 

susceptibles d'etre homologues de la region cible, Le radiologue aura 
alors la possibilite de selectionner manuellement celle qu'il considere 
comme etant effectivement la region homologue, 

Ainsi qu'on Ta 6voque ci-avant, on peut utiliser dans le traitement 

10 d'appariement, un traitement de comparaison specifique au lieu d'un 

traitement de correlation. 

Ainsi, comme illustre plus precisement sur la figure 4, on 
determine (etape 409) pour la region cible contenue a Tinterieur de la 
fenetre cible 6 Id, et a partir des niveaux de gris des pixels contenus 

15 dans cette fenetre cible, une liste de caracteristiques permettant de 

caracteriser en general une region d'interSt. Ces caracteristiques 
peuvent comporter par exemple une valeur representative de 
rSlongation de la region d'interet, une valeur representative de sa 
largeur moyenne, une largeur representative de son niveau de gris 

20 moyen, ou de son gradient moyen d'intensite. 

On effectue la meme operation (etape 410) pour les regions 
d'interets candidates contenues dans des deuxiemes fenetres 61g 
definies a partir de pixels selectionnes dans la zone de recherche. A 
cet 6gard, afin de limiter le temps de traitement, les pixels 

25 selectionnes pour la d6termination des deuxiemes fenetres englobant 

les regions candidates, sont choisis parmi les pixels de la zone de 
recherche presentant des maxima de niveaux de gris dont la dynamique 
est superieure a un seuil predetermine. 

A partir de ces differentes listes de caracteristiques ainsi 

30 obtenues, on effectue un traitement de comparaison 411 entre la liste 

de caracteristiques affectees k la region cible et chacune des listes de 
caract6ristiques affect6es a chacune des regions candidates. 

En pratique, bien qu'on puisse envisager un autre type de 
traitement de comparaison, on utilisera de pr6ference un traitement de 

35 comparaison prevoyant la minimisation d*une distance, par exemple la 
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distance eudidienne, entre la liste des caracteristiques de la region 
cible, consid6ree alors comme un vecteur, et une liste de 
caracteristiques d'une region candidate, consideree egalement comme 
un vecteur. 

5 On pourrait egalement prevoir, a titre d'exemple, de determiner la 

valeur moyenne des valeurs absolues des differences entre les valeurs 
homologues de deux listes de caracteristiques, et de rechercher celle 
des valeurs moyennes qui est la plus faible pour en deduire la region 
d'interSt homologue. 

10 La selection de la region homologue prevue dans l'6tape 405 

s'obtient done par cette variante en minimisant une distance entre des 
listes de caracteristiques plutot qu'en maximisant une correlation 
normee ou une difference normee. 

L'invention prevoit encore certaines ameliorations. 

15 Ainsi, dans la variante de I'appariement prevoyant un traitement 

de correlation, on peut limiter la zone de recherche aux pixels de la 
zone epipolaire presentant des maxima de niveaux de gris dont la 
dynamique est superieure a un seuil predetermine. 

En outre, il est preferable, notamment pour s*affranchir des 

20 problemes de fond d'image, d'effectuer un filtrage prealable des 

images stereotaxiques du type dit "chapeau haut de forme" (ou 
transformation top-hat en langue anglaise). Ce genre de filtrage est 
connu de Thomme du metier et celui-ci pourra se referer pour plus de 
details a Tarticle de J. SERRA Image Analysis and Mathematical 

25 Morphology, Vol. 2, Academic Press 1988. 

Ce genre de filtrage permet d'extraire d*une image numerique les 
zones Claires et etroites et permet de s*affranchir du fond. 

Les etapes du procede selon Tinvention qui vont maintenant etre 
decrites, et prevoyant notamment une reprojection sur Timage de 

30 centrage, s*apparentent ainsi notanmient conune etant une verification 

des resultats de Fappariement precedent entre I'image cible et la 
deuxieme image, et permettent ainsi de minimiser encore les risques 
d^erreurs d'appariement. 

Ceci etant, bien que la suite du proc6de qui va maintenant Stre 

35 decrite en se referant plus particulierement aux figures 4 et 6, prevoit 
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la possibilite d'utiliser des traitements de correlation ou de 
comparaison du type de ceux qui viennent d'etre decrits, ce procede 
selon rinvention utilisant une reprojection sur Timage de centrage, est 
independant de rappariement prealable utilise entre Timage cible et la 
deuxieme image. 

On suppose maintenant (figure 6) que le procede selon Tinvention 
a permis de selectionner sur Timage de gauche 12 deux regions 
d'interets candidates 5g et 20g susceptibles d'etre la region homologue 
de la region cible 5d apparaissant sur Tinnage cible 10. 

A partir du pixel cible de la region dlnt6ret cible et du pixel 
candidat ayant donne lieu a la selection de la r6gion d'interet 
candidate 5g, on dfitermine (etape 406), par un calcul trigonom6trique 
de stereotaxic, et connaissant la geometrie du dispositif de stereotaxic, 
les coordonnees tridimensionnelles du point correspondant de Tobjet 
d'interet candidat correspondant a ces deux regions 5d et 5g. Puis, par 
un calcul trigonometrique de stereotaxic inverse, on determine les 
coordonnees dans Timage de centrage du pixel projete correspondant k 
la projection du pixel de Telement candidat dans cette image de 
centrage. 

On remarque done que, dans le cas present, puisque la region 
candidate 5g est effectivement la region homologue de la region 5d, la 
position du pixel projete dans Timage de centrage 1 1 correspond 
effectivement a une region d'interet projete 5p qui est effectivement la 
projection dans cette image de centrage de relement d'interet 5. 

Par contre, par ces memes calculs, la region d'interet candidate 
20g va conduire a Tobtention d'une position projetee 20p dans Timage 
de centrage qui, dans le cas present, ne se superposera a aucune region 
d'interet physiquement presente dans cette image de centrage. 

Aussi, un traitement de correlation entre la fenetre cible contenant 
la region cible 5d et un voisinage du pixel projete contenant la region 
d'int6ret projetee 5p va donner une valeur de correlation sup6rieure h 
la valeur de correlation obtenue entre la fendtre cible 5d et un 
voisinage du pixel projete 20p. II en aurait 6t6 de meme, meme si le 
voisinage du point projete 20p avait correspondu k une micro- 
calcification projetee autre que celle correspondant effectivement a la 
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microcalcification 5. 

Ce precede permet done de selectionner parmi les deux candidats 
5g et 20g la region 5g qui est effectivement la region homologue de la 
region 5d. 

5 Bien entendu, la conclusion aurait 6te la mSme en utilisant un 

traitement de comparaison 412 au lieu d'un traitement de correlation 
407. 

L'utilisation d'un voisinage du pixel projete, typiquement 
quelques pixels autour du point projete, permet de prendre en compte 

10 les erreurs eventuelles dues a un bouge de Tobjet entre Texposition a 

0°et les expositions ±15°. 

D'une fa9on plus generale, on peut ordonner Tensemble des 
regions candidates obtenues dans la deuxifeme image dans Tordre 
decroissant de leur probabilite a etre la region homologue de la region 

15 cible. Ainsi, la region ay ant par exemple donne le maximum de 

correlation avec la dynamique la plus forte sera affectee du rang 1 et 
ainsi de suite par ordre decroissant. On peut effectuer ce meme 
ordonnancement pour I'ensemble des regions projetees obtenues sur 
I'image de centrage. On selectionnera alors la region homologue de la 

20 region cible comme etant celle dont le produit des rangs respectifs 

dans les deux ordonnancements est minimum. 
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REVENDICATTONS 

1. Proc6de de localisation d'un element d'interet contenu dans un 
objet tridimensionnel a partir des positions de regions d'interets 
homologues correspondant audit element d'interSt et apparaissant dans 
un ensemble d'images stereotaxiques dudit objet, comportant une 

5 selection dans une premiere image stereotaxique (10) d'une premiere 

region d'interet (5d) dite "cible", ainsi qu*un appariement de ladite 
premiere region avec une deuxieme region d'interet (5g) homologue de 
la premiere et apparaissant dans une deuxieme image st6reotaxique 
(12), caracterise par le fait que, les images stereotaxiques etant 

10 num^risees, on selectionne (400) dans Tetape de selection, un" pixel 

cible (60d) dans ladite region d'interet cible (5d), et dans Fetape 
d*appariement, on genere (401), autour du pixel cible selectionne, une 
fenetre cible (6 Id) de caracteristiques dimensionnelles choisies et 
contenant ladite region d'int^rSt cible, on determine (402) dans la 

15 deuxieme image, selon un critfere de selection predetermine, un 

ensemble de pixels, on genere (403) autour de chaque pixel 
selectionne une deuxieme fenetre de memes caracteristiques 
dimensionnelles que ladite fenetre cible, on effectue un traitement de 
correlation (404) entre les niveaux de gris des pixels de chaque 

20 deuxieme fenetre et les niveaux de gris des pixels de la fenetre cible, 

de fa9on a obtenir une valeur de correlation pour chaque deuxieme 
fenetre, et on identifie la region d'interet homologue de la region 
d'intorSt cible h peirtir de Tanalyse de Tensemble des valeurs de 
correlation ainsi obtenues, de fagon a minimiser les risques d'erreur 

25 d'appariement entre les regions d'interet homologues. 

2. Procede selon la revendication 1, caracterise par le fait que 
Tanalyse des valeurs de correlation obtenues comporte la selection 
d'un nombre choisi de maxima ou de minima de corr61ation, la 
selection de la region d'int6ret homologue s'effectuant parmi celles 

30 dont la valeur de corr61ation associee est I'un de ces maxima ou 

minima de correlation. 

3. Proc6de selon la revendication 2, caracterise par le fait que 
I'analyse des valeurs de correlation obtenues comporte en outre une 
determination de la dynamique des maxima ou minima selectionn6s, et 
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la comparaison des valeurs de dynamique obtenues a un seuil. 

4. Precede selon Tune des revendications 1 a 3, caracterise par le 
fait que le traitement de correlation comporte un traitement de 
correlation normee. 

5 5. Procede selon Tune des revendications 1 a 3, caracterise par le 

fait que le traitement de correlation comporte un traitement de 
difference normee. 

6. Procede de localisation d*un element dlnteret contenu dans un 
objet tridimensionnel a partir des positions de regions d'interSts 

10 homologues correspondant audit element d'interet et apparaissant dans 

un ensemble d'images stereotaxiques dudit objet, comportant une 
selection dans une premiere image stereotaxique (10) d'une premiere 
region d'interet (5d) dite "cible", ainsi qu'un appariement de ladite 
premiere region avec une deuxieme region dlnteret (5g) homologue de 

15 la premiere et apparaissant dans une deuxieme image stereotaxique 

(12), caracterise par le fait que, les images stereotaxiques 6tant 
numerisees, on selectionne(400) dans I'etape de s61ection, un pixel 
cible dans ladite region d'interet cible, et dans Tetape d'appariement, 
on genere (401), autour du pixel cible selectionne, une fenetre cible 

20 (6 Id) de caracteristiques dimensionnelles choisies et contenant ladite 

region d'interet cible, on determine (409), pour ladite region d'interet 
cible, a partir des valeurs de niveaux de gris des pixels de ladite 
fenetre cible, un premier ensemble de valeurs numeriques 
respectivement representatives de caracteristiques predetemninees 

25 caracterisant toute region d'intergt, et on stocke ce premier ensemble 

de valeurs numeriques, on determine (402) dans la deuxieme image, 
selon un critere de selection predetermine, un ensemble de pixels 
appartenant respectivement a des regions d'interet dites "candidates" 
apparaissant dans la deuxieme image, et on genere (403) autour de 

30 chacun de ces pixels selectionnes une deuxieme fenetre contenant 

ladite region d'interet candidate correspondante, on determine (410), 
pour chaque region d'interet candidate, a partir des valeurs de niveaux 
de gris des pixels de la deuxieme fenetre associee, un deuxieme 
ensemble de valeurs numeriques respectivement representatives 

35 desdites caracteristiques pred6terminees, on stocke ce deuxieme 
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ensemble de valeurs numeriques, et on identifie la region d'int^ret 
homologue de la region d'interet cible h partir d'un traitement de 
comparaison (411) entre ledit premier ensemble de valeurs numeriques 
et chaque deuxieme ensemble de valeurs numeriques, de fafon h 
minimiser les risques d*erreur d*appariement entre les regions d'interet 
homologues. 

7. Procede selon la revendication 6, caracterise par le fait que 
lesdites caracteristiques predetenninees comprennent des 
caracteristiques de forme, de contraste, de gradient. 

8. Procede selon la revendication 6 ou 7, caracterise par le fait 
que le traitement de comparaison comporte un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

9. Procede selon Tune des revendications precedentes, caracterise 
par le fait qu'on determine dans la deuxieme image stereotaxique une 
zone epipolaire contenant au moins le segment epipolaire relatif au 
pixel cible, et par le fait que les pixels de la zone epipolaire 
contiennent les pixels selectionnes. 

10. Procede selon Tune des revendications precedentes, 
caracterise par le fait que lesdits pixels selectionnes sont choisis 
parmi les pixels de la deuxieme image presentant des maxima ou 
minima de niveaux de gris dont la dynamique est superieure a un seuil 
predetermine. 

11. Precede selon rune des revendications prcCcdentes, 
caracterise par le fait que Tappariement comporte prealablement un 
filtrage des images stereotaxiques, en particulier du type dit "chapeau 
haut de forme". 

12. Procede de localisation d'un 61ement d'interet contenu dans un 
objet tridimensionnel h partir des positions de regions d'int6rets 
homologues correspondant audit element d'interet et apparedssant dans 
un ensemble d'images ster6otaxiques dudit objet, caract6ris6 par le fait 
que, les images stereotaxiques 6tant num6risees, il comprend une 
selection dans une premi&re image st6r6otaxique (10) d'une premiere 
region d'int6ret (5d) dite "cible", une s61ection (405), dans une 
deuxieme image st6reotaxique (12), k psirtir d'un premier appariement 
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automatique, d'au moins une deuxieme region d*interet (5g) suceptible 
d'etre homologue de la region cible, la determination de la position 
spatiale d'un pixel candidat d'un elennent d'interet candidat 
correspondant k ces deux regions dlnteret, la determination (406) de 

5 la position dans une troisieme image stereotaxique, d'un pixel projete 

correspondant a la projection dans cette troisieme image stereotaxique 
du pixel candidat, et un deuxieme appariement automatique entre la 
region d'interSt cible et un voisinage du pixel projete definissant une 
region d'interet projetee, de fa9on a minimiser les risques d'erreur 

10 d'appariement entre les regions d'interet homologues, 

13. Procede selon la revendication 12, caracterise par le fait qu'on 
selectionne dans la deuxieme image stereotaxique, a partir dudit 
premier appariement automatique, un ensemble de deuxiemes regions 
susceptibles d'etre homologues de la region cible, on determine la 

15 position spatiale de chaque pixel candidat correspondant a chaque 

paire de regions d'interet form^e par la region cible et Tune des 
deuxiemes regions, ainsi que la position dans la troisieme image de 
chaque pixel projete correspondant, et on effectue le deuxieme 
appariement automatique entre la region cible et un voisinage de 

20 chaque pixel projet6. 

14. Procede selon la revendication 13, caracterise par le fait qu'on 
ordonne Tensemble des deuxiemes regions obtenues dans I'ordre 
decroissant de leur probabilite a etre la region homologue de la region 
cible, on ordonne I'ensemble des regions projetees obtenues dans 

25 I'ordre decroissant de leur probabilite h Stre la region homologue de la 

region cible, et on selectionne la region homologue de la region cible 
comme etant celle dont le produit des rangs respectifs dans les deux 
ordonnancements est minimum. 

15. Procede selon Tune des revendications 12 ^ 14, caracterise par 
30 le fait que le premier appariement automatique est un appariement tel 

que defini dans Tune des revendications 1 a 11. 

16. Procede selon Tune des revendications 12 a 15, caracterise par 
le fait que le deuxieme appariement automatique comporte un 
traitement de correlation tel que defini dans Tune des revendications 1 

35 as, entre le voisinage de chaque pixel projete et une fenetre cible 
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contenant la region d*int6ret cible, ou un traitement de comparaison, 
tel que d6fini dans Tune des revendications 6 a 11 entre la region 
d'interet cible et chaque region d'interet projetee. 
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